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Abstract 

The strong coupling half-filled Kondo lattice model is an important example of a strongly in- 
teracting dense Fermi system for which conventional Fermi gas analysis has thus far failed. We 
remedy this by deriving an exact transformation that maps the model to a dilute gas of weakly 
interacting electron and hole quasiparticles that can then be analyzed by conventional dilute Fermi 
gas methods. The quasiparticle vacuum is a singlet Mott insulator for which the quasiparticle 
dynamics are simple. Since the transformation is exact, the electron spectral weight sum rules are 
obeyed exactly. Subtleties in understanding the behavior of electrons in the singlet Mott insula- 
tor can be reduced to a fairly complicated but precise relation between quasiparticles and bare 
electrons. The theory of free quasiparticles can be interpreted as an exactly solvable model for a 
singlet Mott insulator, providing an exact model in which to explore the strong coupling regime of 
a singlet Kondo insulator. 
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Introduction. A Fermi liquid is described by starting with an ideal non-interacting 
Fermi gas. We then adiabatically "switch on" the interactions between particles. States of 
the interacting system are then identified with those of the non-interacting system. Fermi 
quasiparticles have the same relation to Fermi liquid ground state as the electron operators 
do to the vacuum; if | ty G ) is the Fermi liquid ground state, there should exist a quasiparticle 
creation operator cjg so that = cp \ ) , { c a,c^} = 5 a p and Hc\\ ty G ) = e a c a\ )■ 
The spectral weight sum-rule follows directly from this. 

Although the spectral weight sum rule is so fundamental it has been surprisingly difficult 
to implement in approximate theories of the Kondo lattice model, jl| The reason is that 
the ground state is a dense system of strongly interacting Fermions for which it has been 
difficult to obtain results using traditional weak coupling methods. This work demonstrates 
an exact transformation that can be utilized to map the dense strongly interacting singlet 
Mott insulator to a weak coupling dilute Fermi gas. This allows conventional Fermi gas 
techniques to be applied to this class of models. 

We begin by considering a lattice of sites each consisting of conduction band "c" electrons 
that can hop between neighboring sites and core "f" electrons that are confined. At half 
filling a total of two electrons tend to be located on the site. A potential favors single 
occupancy of the "f" electron, forcing the second electron into the conduction band. The 
extended model consists of a lattice of such sites with only the "c" electrons hopping between 
neighbors. 

There are four states with one "c" and one "f" electron present. The resulting triplet- 
This behavior is most easily encoded in the Kondo lattice Hamiltonian |3[ 

h klm = t ErA4A r ) c cA r ') + cc)+ . (i) 

E r (JS c (r)-S f (r) + U f (n f (r)-lf) 

The atomic ground state is given by | ^/ G ) = ^ (4 t c /i ~ c c .1 4?) I ). The full ground 
state is a product of such two-electron singlets at each site and is nondegenerate with finite 
gap to the excited states. The gap will persist for small values of the hopping t and the local 
conduction electron will then simply make virtual excursions to the neighboring sites. The 
system will remain a singlet Mott insulator. 

Before describing quasiparticles for nonzero t we begin with the limit of zero hopping. For 
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the Hubbard impurity, there are four eharged spin half exeitations on a single site, consisting 
of each of the four states with three electrons. |3j We first attempt to define a quasiparticle 
operator simply by which one might hope adds an extra conduction electron to the 
singlet. 

This ansatz immediately leads to problems. First we find c\ * | ) = 
~ 75 c c T c c | C /T I ^ )■ Although the operator indeed creates a state with an extra charge, 
this state is not normalized. 

A second problem is the following inequality: c c ^ | ) = j_ | ) ^ which shows 
that the annihilation operator does not annihilate the atomic ground state. Finally, we note 
that c c | | ^>q ) = — ^75 c c j, c / t I ) 7^ I ) so that this simple annihilation operator 
does not reconstruct the original singlet ground state. We conclude that none of the condi- 
tions we require of a proper quasiparticle operator is fulfilled with this ansatz. This simple 
analysis reveals the source of the difficulty; the ground state is "entangled" in the sense that 
it cannot be created from the true vacuum by a simple product of bare electron operators. 

Constructing correct local quasiparticles. 

We now discuss how to define proper quasiparticle operators. The local Hilbert space 
consists of sixteen states which we label from 1 to 16. First we list the eight states with 
even particle number ordered first by increasing particle number. Degeneracies are split by 
total spin, and remaining degeneracies are split by s z . States with odd particle number are 
similarly assigned labels 9 to 16. 

This "canonical" ordering is shown on an energy versus particle number diagram for 
the Hamiltonian h canonica i = n tot — | (rif — n c ) in Figure On the horizontal axis is 
plotted ntot- On the vertical axis is the energy state labelled with the state of minimum 
z-component of spin for each degenerate multiplet. In Figure [IB is a hypothetical diagram 
with an interacting Hamiltonian h atom i C favoring the Kondo singlet ground state. To avoid 
cluttering the figure, triplet states are not drawn. 

We denote \n) a =\n ) atomic states in Fig. Hp and | n ) c = | n ) canonical states states in 
Fig. [TK. A precise description of these states is needed to continue. 

The singlets are given by | 1 ) c = 1| ), | 2 ) c = 4 T J }, | 3 } c = ^= (4,t 4,1 ~ 
4^ c| iT )| ), | 4 ) c = 4,t c c,||0)> | 5 ) c = 4,t4,| c /,t c /,|l )■ The triplets numbered 
| 6 ) c , | 7 ) c and | 8 ) c are not shown in the diagram. 
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Energy vs particle number 




FIG. 1: In (A), the energy E = ntot — |(^/ — n c ) is plotted vs particle number and labeled with 
the state that has the smallest value of S 2 and S z . A dotted arrow indicates a c-type creation 
operator is used to create the state at the tip of the arrow from the state at the base. A a solid 
arrow indicates an f-type creation operator is used. A solid-dotted combination indicates a singlet 
combination is used. In (B) a similar construction is used in an interacting model. The dotted 
arrow here indicates a hole operator. 

The four states of total spin — | are given by: | 9 ) c = 4 ,| ), | 10 ) c = c\, \ ), | 11 ) c = 
c c,J. c /,T c /,|l ), | 12 ) c = c I,t c c,| c /,jJ ). The states of total spin ~: are | 13 ) c = 
C y ), | 14 ) c = 4 jT | ), | 15 ) c = < t 4 jT 4J ), | 16 ) c = 4 )T 4 4 4 iT | >. The states 
with s z = — | are not labelled in the figure. 

A local spin-symmetric Hamiltonian H is block diagonal in the blocks enclosed by paren- 
theses: ([1], [2,3,4], [5]) ([6]), ([7]), ([8]) ([9, 10][11, 12]) ([13, 14][15, 16]). If particle number is 
also conserved, H is further block diagonal in the subblocks in square brackets. 

We now draw in Fig. [JJ3 a similar diagram for a more complicated Kondo-Lattice like 
Hamiltonian with identical but permuted eigenvalues. We assume it has the Kondo singlet 
ground state and lowest energy charge 2 ± 1 states consisting of zero or two c-electrons 
together with a single "f" electron. Eigenstates in Fig. [TJ3 are labelled according to the 
same ordering scheme. 

Our goal is to construct a unitary transformation U that maps the Kondo singlet ground 
state in Fig. [TB to the vacuum in Fig. [TJA. and at the same time define quasiparticle operators 
that create the states of charge 2 ± 1 from the Kondo singlet. The arrows in Fig. [T]B illustrate 
the action of the quasiparticle operators. A dotted hole creation operator creates the singly 
charged states and or a solid electron operator creates the triply charged states. We thus 
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map the ground state | 3 ) a to the vacuum | 1 ) c = | G ) c i.e. U\ 3) a = | 1 ) c . The low 
energy state of one extra spin down electron | 12 } o must map to | 9 ) c ; U\ 12 ) o = | 9 ) c . 
Similarly we demand U\ 9 ) = | 10 ) . 

Symmetries can now be used to completely specify U. We define Q to be the combined 
particle-hole spin flip transformation. This permutation is an inversion of the picture in Fig. 
[TJ\ through the state | 3 ) c . Let be the 16 x 16 matrix that has all entries zero except 
Gij = 1. It can be verified that Q = Gi j5 + £3,3 + G 5:1 + G 6j6 + G 7j7 + G 8 ,s + £9,12 + Gi ,ii + 
Grs,w + Gu,i5 ~ G 16A3 - G 15>u - G 12 , 9 - G u ,io - G 4:2 ~ G 2A . Signs are chosen to preserve 
spin. 

We next define K to be the transformation that exchanges the two species of Fermions. 
This corresponds to the permutation obtained by exchanging the high and low energy state 
for each value of n tot in Fig. [T]A.. Minus signs are inserted to make the transformation an 
element of the group of continuous rotations. We find K = G\^ + G 2 ,4 + G^ 2 + G 5)5 + G 6j6 + 

G?j + Gs,8 + Cio,9 + ^11^2 + Gi4 5 i3 + Gl5,l6 ~ G*i6,15 _ ^13^4 — G*12,11 ~ Gg^o — (^3,3- 

In order for U to preserve spin and also generate particle and hole operators from the 
original Fermions, we demand U obeys the following: [U, S] = 0, ntot U = U (2 + n e — nh) 
and UK = QU . It can be verified that these constraints together with the demand that U 
be real and unitary results in 

U= G M + G 2 , 3 /v / 2 + -G' 2 ,5/v / 2 + G'3,i +. 
G^/V2 + G^ 5 /V2 + G 5 , 2 + G 6 

,6 + 

Gfj + Gsfi + Gg } io H Gio,l2 + 

^11,11 + Ci2,9 + ^13,14 H £"14,16 + 

^15, 15 + ^16,13 

This defines a new Hamiltonian h' atomic = U}h canonica iU r which has the energy diagram shown 
in Fig [TJ3 . The demand that TP n to t U = 2 + n e — shows that U cannot be continued to 
the identity transformation which is presumably relevant for J = 0. This observation adds 
support to the existence of a quantum phase transition for sufficiently large values of tj J, a 
subject which is outside the scope of this investigation. 

From unitarity follows that c£,o-( r ) = Tj\(? c ^{r)TJ r and & hrT {r) = Uj.c^ a (r)U r preserve the 
Fermi anticommutation relations. Inverting this formula and using the fact that U r \ ) a = 
U r \ 3 ) Q = | 1 } c = | Vac ) c we find that c cf a {r)\ Vac ) c = is equivalent to c e/l(7 (r)| ^f G } a = 
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0. We have therefore obtained a local unitary transformation that maps the original "c" 
and "f" bare electron operators to quasiparticle operators that annihilate the Kondo singlet. 
By construction the operators $ eh {r) either add or remove exactly one charge in each of the 
states and c\ a {r) is identified as an electron quasiparticle operator and & ha (r) as a hole. 



Using the method in Ref. [4] and Ref. 



5| we can explicitly write down the original 



operators c\j- and c c f as a polynomial of Fermion operators c[ h a (r) , c eh a (r) . The result is 
c hd (ri n hA n e>T + t 2 n eA + t x n M n eA + r x n e>T n eA + \ c{ A c{ A - t x n eA - ^ n M ) + 
where r 2 = (1 - 1/^2) and n = |(^2 - 1). For the "f" electrons I find 



c 
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Similar formulas occur for the other spin component. We note that the "f" electron operators 
is constructed by third and higher order electron and hole operators in contrast to the 
conduction electron operators that have a linear coupling to quasiparticles. Deriving or 
even verifying these formulas is formidable without the help of a computer. [s| 

Effective model for quasiparticles. The transformation U is now used to define a 
global transformation by U g i b a i = Yl r Ur- Since U r does not mix states of even and odd 
particle number we obtain {c ia (r), ct (r')} = Sij5 r y verifying that U is indeed a global 
canonical transformation. 

We now express the Kondo lattice model given by Eq. [T] in terms of electron and hole 
quasiparticle operators by replacing the Fermi operators c c / j(T (r) by their representation as 
polynomials in c\ ha (r). 

By construction, the Kondo singlet is the quasiparticle vacuum in the limit t — > 0. For 
small values of t, quasiparticles can hop between adjacent sites. The resultant low energy 
charged charged states that appear in the interacting ground state are therefore described by 
single quasiparticle operators and the unitary transformation guarantees these local states 
also have the correct energy. 

For small t the states with local charge fluctuations will have a small amplitude in the 
interacting ground; the probability that a site is occupied by more than one quasiparticle 



will be accordingly less. The system should therefore be well approximated by a dilute fermi 
gas. 

We thus expand H KLM in quasiparticle operators. Inserting the quasiparticle represen- 
tation for the original fermions, normal ordering then truncating the Hamiltonian to second 

i j2 l \r\[ tt zj-free . Tjinteractinq .,1 

order we fmd[6|J H KLM = H J KLM + H KLM with 

H f K e L M = \J{n c + n f -l) + 

§tE s (-1)*+' (c{ s (r) c{_ s (r') - c e>s (r') c h ^_ s (r)) + 

ktE s ct tS (r)c e!S (r')-c{ s (r)c h>s (r') + 

\t^ e y)c e ^r)-c\y)%sir) + 

- Iffi, (-1)1+* (c, jS (r') c e ,_ s (r) - gt f .(r) <„ s (r')) 

Hklm can then be treated exactly by transforming to momentum space and utilizing 
a conventional Bogoliubov transformation. [5] With a chemical potential we compute the 
eigenvalues exactly. We let e k = cos(A; • i). The quasiparticle spectrum is then given by 



E k = ^A k ±te k ±n where A k = ^1 + (ff ) 2 . 

It is reassuring that Uf does not appear to quadratic order; a J which is not small is 
what is crucial in creating the Kondo singlet and low-energy states. For completeness, the 
entire onsite part of rl KLM is given by n KLM y = 

{n e (n e - l)n h + n e n h (n h -1))(±J+=±U) + . 

(n e (n e -l)+n h (n h -l))(=fJ+lU) + 

{{S e -S h ){\J-U)+n e n h {=^J+\U)) 
High order terms in the hopping become too lengthy to write out here, but are unimportant 
for the small t physics since not only is the prefactor t small, but also, as we shall see , the 
fourth order Fermion operators become unimportant since the Fermi gas is dilute. 

Spectral weights. Since the local electron operator is represented in terms of local 
quasiparticles the spectral weights can be computed for H^hM- an example, the results 
are plotted in in Fig. [2] where the spectral weight for the c-electron A c (k,u) is plotted as 
a function of u several values of k from to tt. The quasiparticle coherent peak is seen 
as a broadened 5 function. The incoherent structure is comprised of three-quasiparticle 
contributions. Not shown are the five-quasiparticle contributions that represent less than 1 
% of the total spectral weight. We observe that A c {k + tt,u) = A c (k, —uj) a consequence of 
the particle- hole symmetry of H KLM . 
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FIG. 2: Square root of the spectral weight for c-electrons A c (k,uj) is plotted for Hj^^ e M for 
t/J = 0.4 and t = 1 in one dimension. To enable plotting the coherent and incoherent parts 
together, A c {k,oS) has been smoothed by a Gaussian in u before taking the square root. Plots are 
shown for k = to k = it. 

Numerical results. The quality of the wave function as an approximation to the Kondo 
Lattice ground state can be measured by the overlap of the ground state with the Hilbert 



space of H KLM . Let a = \ fff Td {l - A k l ) w f A 2 + 0(A 4 ) where A = t/J and d is 



the dimension. The deviation from unity is measured by 5p = ( (1 — rif) 2 ) which given 
exactly by Sp = a 2 (3 — 2a) demonstrating convergence like (t/J) 4 " for small t . We also find 
( n e + rih ) = 4 a which verifies that the number of quasiparticles remains surprisingly small 
even as the Kondo lattice model approaches a phase transition. 

A more challenging quantity to compute is the expectation value of the complete full 
Kondo Lattice Hamiltonian with Uf = in the ground state of H^^ e M . This quantity is 
important since it is a rigorous upper bound to the true Kondo lattice ground state energy. 



( H KLM ) free = =f + 3a + =f a 2 + fa 3 + =f (3 A 2 + 
a 3 /3A 2 (72 + -48v / 2) + a/3 3 rf" 2 A 4 (f + ^ v / 2) + 
a/3A 2 (16 + -8v / 2)+/3 5 rf- 4 A 6 (^ + ^ v / 2) + 

^3 d -2 A 4 ( ^4 + |^2 ) + a 2 /5 3 d -2 A 4 (z |4 + 128^2 ) + 

a 4 /3 A 2 (-36 + 24 y/2) + a 2 (3 A 2 (-52 + 32 y/2) 



We define 0= JJJ Td e 2 /A k (g ) 
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This can be compared to high order series expansions. 7] We find our result is at most 
about 3 % above the series result, converging very rapidly for small values of A. Errors are 
bounded by .61 A 4 for d = 1, 2.25 A 4 for d = 2 and 4.9 A 4 for d = 3 in the entire domain of 
convergence of the power series. 

Conclusions. An exact canonical transformation is derived that maps bare electron 
operators to quasiparticle electron and hole operators. This allows us to approximate the 
Kondo lattice model as a dilute Fermi gas for large values of the Kondo coupling. Properties 
of the Fermi gas can then be computed exactly. The overlap per site between the free Fermi 
gas ground state and the Hilbert space of the Kondo Lattice model is (1 — 0(t/J) 4 ) and 
ground state energies are also accurate to this order. Spectral weights can be computed that 
exactly obey the sum rules. 

Together with the exact relation between quasiparticle operators, the free quasiparticle 
Hamiltonian represents an exactly solvable model of a singlet Mott insulator. The difficulty 
of understanding electrons in this Mott insulator is encoded in the complicated relation 
between electrons and quasiparticles rather than any subtlety in the quasiparticle dynamics. 
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